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for the B — > decay. The latter mixing is calculated here for the first time. We discuss 
general sum rules which have to be satisfied in any model of this type. We emphasise that 
the neutral gauge bosons in question could also significantly affect other fermion radiative 
decays as well as non-leptonic two-body B decays, e'/e, anomalous (g — 2) M and electric 
dipole moments. 
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1 Introduction 

Extensions of the Standard Model (SM) with additional gauge symmetries are of particular 
interest in view of current direct and indirect searches for physics beyond the SM. A feature 
of such theories is the existence of new gauge bosons. These could provide clear deviations 
from the SM predictions and have been studied as possible early discoveries at the LHC. 

Constructions with a single additional U(l) factor correspond to an extra neutral gauge 
boson Z' [1-4] and have received a lot of attention, since they are naturally present in a large 
variety of models, such as Grand Unified Theories (GUTs), higher-dimensional models, 
little or composite Higgs models, superstring constructions. Hence, the predicted Z' mass 
spans a wide range of values, from the high GUT scale down to the TeV scale, having 
different phenomenological consequences. On the other hand, it is difficult to determine 
model independent bounds on the Z' mass, since its relation to observables depends strongly 



- 1 - 



on the mixing of the Z' to the SM the strength of the new gauge-coupling constant, the 
charges of the SM fermions and Higgs doublets, the presence of additional exotic fermions 
necessary to cancel the anomalies, etc.. 

Still, there exist bounds on the Z' mass and its mixing with the SM Z from precision 
electroweak data, direct and indirect searches at the Tevatron, and interference effects 
at LEP2 in specific Z' models (see [1] and references therein). Assuming an electroweak 
strength for the new gauge coupling constant and family-universal SM fermion charges, 
most of the canonical Z' models have predictions in agreement with the experimental 
constraints, for Z 1 masses not too much below 1 TeV. However, considering leptophobic Z' 
models or models with family non-universal SM fermion charges, the Z' mass is allowed to 
be even smaller. 

Beyond the simple U(l) group, extended gauge symmetries have also attracted a lot of 
attention, especially to explain the observed pattern of fermion masses and mixings (some 
of the earliest papers on this subject are [5-15]). A common feature of these models is the 
appearance of dangerous flavour-changing-neutral-current (FCNC) contributions, mediated 
by the new flavour gauge bosons. Hence, a high new physics (NP) scale, usually larger 
than 1000 TeV, must be adopted to provide the necessary suppressions. This prevents any 
NP direct observation at the present and future colliders. 

An alternative approach to provide an explanation of the pattern of fermion masses is 
considering global flavour symmetries; in particular discrete symmetries have been widely 
studied in the last years (for a review see [16]). It is interesting to notice that the FCNC 
contributions are usually less constraining in these models [17-27] compared to gauge 
flavour models, opening the possibility for a lower NP scale and correspondingly for new 
particles beyond SM that can be detected in present and future colliders. 

Recently, a series of papers [28-31] appeared suggesting the possibility of explaining 
fermion masses and mixings with a NP scale of a few TeV, but with extra gauge symmetries. 
In these models, the gauge flavour symmetry is the product of non-Abelian U(3) factors. 
For each gauge-group generator there is a new gauge boson, with no color or electric charge. 
In addition, new exotic fermions have to be introduced to cancel the anomalies of the new 
gauge sector. 

The analyses in [28-31] and most papers to date studied almost exclusively the impact 
of flavour changing heavy neutral gauge bosons on AF = 2 transitions and rare B and 
K decays, leaving their impact on the radiative decay B —> X s j aside. As the latter 
provides generally a very strong constraint on extensions of the SM, it is important to 
investigate whether the presence of neutral gauge bosons with flavour violating couplings 
and new exotic heavy fermions would put the models in question into difficulties and 
whether generally the role of such new gauge bosons and new fermions in the B — > X s ^f 
decay is significant. 

Having this in mind, we point out in this paper two new contributions to the B — > 
X s j decay, that to our knowledge have not been considered in the literature so far. The 
first one originates from one-loop diagrams mediated by heavy neutral gauge bosons and 
exotic quarks with electric charge —1/3, when flavour-violating left-handed and right-handed 
couplings of these gauge bosons to SM and exotic quarks are present. Then, this new 
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physics contribution to the Wilson coefficient of the dominant operator is enhanced by 
m' F I m& with m& being the mass of the 6-quark and m' F the mass of the exotic fermion, in 
general significantly heavier than SM ones. 

Analogous enhancement is known from the study of left-right symmetric models [32— 
34], but the absence of new heavy fermions and the smallness of Wl — Wr mixing makes 
this contribution, mediated dominantly by charged SM gauge bosons with right-handed 
couplings, sufficiently small to agree with data. 

Within the littlest Higgs model with T-parity, similar contributions to the B — > X s ^f 
decay, mediated by heavy neutral gauge bosons and heavy mirror fermions, have already 
been considered [35]. However, in this context only left-handed currents are present and 
in addition the diagrams in question are strongly suppressed by the mixing angles in the 
T-odd sector and by the GIM mechanism, so that these contributions turn out to be very 
small. 

The second new contribution stems solely from the heavy neutral gauge bosons, inde- 
pendently of the presence of exotic fermions in the theory. These gauge bosons are inte- 
grated out at the scale of the order of their masses, in general much higher than /xj, « m^. 
The resulting neutral current- current dimension-six operators affect the B — > X s ^y rate 
through QCD mixing with the magnetic dipole operator responsible for the B — > X s ^f de- 
cay. This contribution is present when the new neutral gauge bosons have direct couplings 
to the b and s quarks. 

As we shall demonstrate explicitly, these new contributions imply constraints on the 
parameters of models, containing neutral flavour gauge bosons with couplings described 
above and with masses within the reach of the LHC. We shall also identify the condi- 
tions under which the first of these contributions turns out to be negligible, such as very 
small mixings among SM and exotic quarks, naturally following from a fermion see-saw 
mechanism. 

We point out that similar contributions could be considered in fj — > f^, like \i — >• 
and t — > cj, where / is any SM fermion, mediated by exotic fermions with the same electric 
charge as /. The same comments apply to flavour conserving observables like (g — 2)^ and 
EDMs, in which dipole operators play the dominant role. Moreover, our QCD analysis of 
the mixing among neutral current-current and QCD-penguin operators is also applicable 
to other processes, such as non-leptonic two-body B decays, and other observables, like 
e'/e. 

We describe the general framework and fix the notation that will be used throughout 
the rest of the paper in Sect. 2. In Sect. 3 we present the analytical results for the exotic- 
fermion contribution to the B — > X s j decay in the context of only one neutral flavour 
gauge boson and one exotic quark with electric charge —1/3. Subsequently, we generalise 
our results to an arbitrary number of such gauge bosons and quarks. 

In Sect. 4 we discuss the second new contribution to the B — > X s j decay. We extend the 
SM QCD renormalisation group (RG) analysis by taking into account all the new neutral 
current-current operators, generated through the exchange of a single neutral gauge boson 
Ajj- To this end we calculate the relevant anomalous dimension matrices at the leading 
order (LO). 
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In Sect. 5 we present our model-independent results in a general form, that can be 
easily applied to all models with flavour-violating neutral gauge bosons, like Z'-models or 
models with non-Abelian gauge flavour symmetries. We identify three major classes of 
such models; for each we cast our results into simple formulae, that depend only on the 
masses and couplings of the new particles. Finally, we illustrate how these formulae can be 
translated into constraints on the parameter space of representative toy-models. In Sect. 6 
we summarise and conclude. 

2 Notation and fermion see-saw mechanism 

We focus on models in which the new physics contributions to5-> X s ^f discussed in the 
previous section are present. In what follows, we describe a general framework and we fix 
the notation. 

We consider the general case in which a product of k > 1 gauge flavour groups, 
Qf = Qx x . . . x Qfr, is added to the SM gauge group. Depending on the dimension di of 
each symmetry factor, d\ + . . . + d^ new flavour gauge bosons enrich the particle spectrum. 
In general, the SM fermions ip transform under Qf and consequently anomalies from gauge 
or mixed gauge-gravitational triangles can arise. In this case, it is possible to introduce a 
new set of exotic fermions suitably transforming under the SM and gauge flavour groups 
to obtain an anomaly-free theory. 

On the other hand, the presence of these new exotic fermions is a quite general feature 
of flavour models, when considering renormalisable ultraviolet completions. This holds not 
only for gauge flavour symmetries, that we are discussing in this paper, but also for global 
flavour symmetries. Indeed, in these models the new exotic fermions and their mixing with 
the SM ones account for the mass spectrum of the SM fermions (see [36] and references 
therein). This typically happens through a fermion see-saw mechanism where the small 
masses of the SM fermions are explained by the large masses of their corresponding exotic 
fermions. Without loss of generality, we consider only the case of S'L r (2)i-singlet exotic 
fermions. Hence the see-saw Lagrangian can be schematically written as 

£ see -saw = M i£ R + i> L V R (j) + y L M 2 D 4> R + h.c. , (2.1) 

where ipL,R (^r,l) comprise all left- and right-handed SM (exotic) fermion fields; their 
dimensions also define the dimensions of the matrices M, Yf and M®. In the same 
expression, <f> represents the SM Higgs field. 

This see-saw Lagrangian can be made invariant under Qf if both SM and exotic 
fermions belong to the same representation of Qf and M, Yf, and are proportional 
to the identity matrix in the flavour space. This case, however, is phenomenologically not 
interesting, because then the flavour symmetry is exactly preserved and it is not possible 
to reproduce the observed pattern of SM fermion masses and mixings. 

On the other hand, a more realistic situation appears when we impose the invariance 
under Qf at high scales, by introducing additional degrees of freedom, and then break the 
flavour symmetries at lower scales to provide the correct SM fermion masses and mixings. 
In the see-saw Lagrangian, the invariance under Qf and its breaking mechanism can be 
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achieved in various ways: we can promote M, Y®, and MP , or only some of them, either 
to spurion fields a la MFV [37-40] , or to dynamical scalar fields with non- vanishing vacuum 
expectation values (VEVs) [41], or to fermion condensates [42]. In what follows, we simply 
treat the terms with M, Y\ ■, and MP, or only some of them, as flavour violating terms 
originating from some non-specified dynamics at the high scale [43]. 

In the simplest phenomenologically interesting scenario, we discuss, the fields tpi (V'-r) 
and ^ r {^! l) transform under the same representation of a single flavour group Qf. The 
see-saw Lagrangian is the same as in Eq. (2.1), with the only flavour symmetry breaking 
parameter being M, while and MP are family universal parameters. In general, the 
square matrix M has non-diagonal hierarchical entries in order to provide the correct 
description of the SM fermion masses and mixings. Before finding the mass eigenvalues of 
the SM and exotic fermions, it is useful to diagonalise M. In general, M is diagonalised 
by a bi-unitary transformation 

(2.2) 



M = V} MV R , 



where Vl,r are unitary matrices and M is a diagonal matrix. After the absorption of Vl 
(Vr) in (^ r) and tpR (tPl) by a field redefinition and the breaking of the electroweak 
symmetry, the resulting mass matrix is as follows: 







YP(<f>) x 1 



MP x 1 



M 




(2.3) 



Reabsorbing possible unphysical phases, we can now define the mass eigenstates $l,r 
and Fl^r by performing orthogonal rotations on the flavour eigenstates ipL.R an d *&l,r- 
For each generation we may write 




cos 9 L)R -smO L)R 
sin 6>l,r cos 6 L)R 




(2.4) 



Denoting with m f k the mass of the light fermion and with vn! Fh that of the heavy exotic 
fermion Fk and using for brevity the notation MP = YP ((/)), what follows is a direct inverse 
proportionality between m,f k and m' F , as in the usual see-saw mechanism: 



m h m' F =MPMP 



(2.5) 



We can express such masses in terms of the flavour symmetry breaking terms: 



sin0,R, sm6 L , 



-M k , 



/fc cos 2 9 Rk - sin 2 e Lk 
where a straightforward calculation gives 



m F k 



cos6 Rk cos6 Lk 
cos 2 9 Rk - sin 2 8 Lk 



-Mk 



(2.6) 



sin 9 L , 



ft 
MP 



k MPm' Fk -M»m fl 



III' 



12 



m. 



sin Or, 



MPm' F -Mfm fi 



\MP 



nr 



12 



m 



■ (2-7) 



Ik 
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These results are exact and valid for all the fermion generations. However, taking 
the limit in which m' Fk 3> m/ fe and Mj° ~ = M D , we find simple formulae that 
transparently expose the behaviour of the previous expressions. In this limit we find 



D\2 



(M u ) 



m' Fk « M k 



smU Lk 



sm( 



'R k 



l m fk 



(2. 



as in the usual see-saw scheme in the limit M k 3> M^ 2 - Notice that these simplified 
relations are valid for all the fermions, apart from the top quark for which the condition 
M3 S> Mj° 2 is not satisfied and large corrections to Eq. (2.8) are expected. 

When the flavour symmetry is broken, the flavour gauge bosons develop masses pro- 
portional to the flavour breaking parameters. It is not possible to enter into details without 
specifying a particular model and in the following we simply assume that the neutral gauge 
boson masses are controlled by the parameters as for the exotic fermions. 

The part of the Lagrangian containing the kinetic terms and the usual couplings among 
fermion and gauge bosons is given by 



c kin = i (^ip^ + ^ipy) . 



(2.9) 



Since we are only interested in discussing the couplings between fermions and the new 
gauge bosons, in writing the relevant covariant derivative we do not consider the coupling 
terms with the SM gauge bosons but only the following term 



(2.10) 



In this expression, gn is the flavour-gauge coupling, T-^ are the generators of the flavour 
group Qf in the representation 1Z, and Ajj a the corresponding flavour gauge bosons. The 
resulting Feynman rules for a single gauge boson Ah in the fermion-mass eigenbasis, T = 
{F, /}, are given by the following relations: 



with 



A 



H 




^I^PlCl + PrCr) 



A H fifj : < 



C L 
Cr 



■ cos ti Li 

■ cos 8r 



V r T k V R) .. cos Lj + sin 9 Li {ylT%V L ) sin 9 L 



Vl T k V l) . . cos 6 Rj + sin Or, (vJt^Vr) . sin I 



Ah Fi Fj : < 



'C L = sm9 Li {y ] R T^V R ) i ^9 Lj + C o S 9 Li (v£z%Vb) ms9 Lj 
C R = sin 9 Ri (ylT^V L ) . sin 9 Rj + cos 9 Ri (v ] R T%V R ) . cos 9 Rj 



(2.11) 



(2.12) 



(2.13) 
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A H fiFj : < 



' c L = cos e Li (v ] r t%v r ) sin e Lj - sin e Li (ylT%y L ) . . cos e L 

' C K = cos0 fli (vJ^Vb) sin0 R . - sin^ (v£t£Vr)^ cos 



(2.14) 



f C L = sin £i ( v£z£ Vr) cos l - cos 9 Li (V^ V L ) sin L 

An^i/i ■{ ) i \ ij / t \ j ( 2 - 15 ) 

[ C fl = sin ^ (v£Z£VlJ cos Hj - cos 0^ (V^Vr) sin Hj 

A few comments are in order: 

1) These results are valid for both Abelian and non-Abelian flavour gauge symmetries. 
In particular, barring accidental cancellations, the couplings described in Eqs. (2.14) 
and (2.15) are always non- vanishing. The only exception consists in all those Z' 
models in which the SM fermions have universal charges. In these cases the flavour 
gauge symmetry Qf is a simple U(l) and the charge universality assures that oc 1. 
Hence, the couplings in Eqs. (2.14) and (2.15) identically vanish and the couplings 
in Eqs. (2.12) and (2.13) turn out to be flavour conserving. 

2) When the flavour symmetry is a product of different groups, Qf = Q±x...x Q^, 
and fermions transform non-trivially under different factors of Qf, Eqs. (2.12)— (2.15) 
change. In particular, if tpL and ipR (\J/r) transform under two distinct groups, 
then only one term in the right-hand side of the previous equations is present (for 
example in Eq. (2.12), Cx could have either the term cosO^ (v^T^VrJ cos 6^. or 

sin (9/^ {^L^^yLj s i n ^Lj i Du t n °t both terms simultaneously). As a result, even for 

the simple case Qf = C7(l) 2 with universal fermion charges, the couplings described 
in Eqs. (2.14) and (2.15) are non- vanishing. Notice that in this particular case all 
the couplings in Eqs. (2.12)-(2.15) turn out to be flavour conserving. 

3) When the SM fermion masses originate from the see-saw mechanism with heavy 
exotic fermions, sm0L k and sin0R k are small (cf. Eqs. (2.7) and (2.8)). Hence, the 
couplings described in Eqs. (2.14) and (2.15) are all suppressed. In the particular 
case of m' Fk 3> mj t , this suppression is roughly given by ^Jm f k / m' Fk . 

4) When a theory is anomaly free without the necessity of introducing extra exotic 
fermions, only the expressions in Eq. (2.12) apply. Still the couplings can be flavour 
violating. This is the case of U(l) models in which either the charge assignment of the 
SM fermions is such that no anomalies arise or the Green-Schwarz mechanism [44] is 
implemented into the theory to compensate for the anomalies. In the specific case of 
Z' models with universal fermion [/(l)-charges, the couplings are flavour conserving, 
as already discussed in 1). 

As we will show in Sects. 3 and 4, in virtually any gauge model, in which SM fermion 
masses are explained through a see-saw mechanism with heavy exotic fermions, both new 
contributions we introduced in Sect. 1 are present and have potentially interesting effects. 
Their impact depends on the details of the specific model, according to the above discussion. 
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On the other hand, for completeness we will also consider models with gauge sym- 
metries beyond those of the SM that provide either an explanation of the fermion mass 
patterns through a different mechanism than the see-saw or no explanation at all. In these 
cases, heavy exotic fermions might not be present in the particle spectrum, as in 4), and 
the corresponding contributions are absent. However, when present, the couplings listed 
in Eqs. (2.12)— (2.15) have in general a completely different structure. In particular there 
might be no suppressions of the exotic fermion contributions due to sin#£ fc and sin 6 R k . 

The presence of exotic fermions can modify the couplings of the SM W and Z bosons 
to fermions with respect to their couplings within the SM. As a consequence, non-unitary 
quark and lepton mixing matrices and flavour violating Z couplings are usual results of such 
changes, which are specific to the considered gauge flavour model. Also these modifications 
affect only the Wilson coefficients and do not generate further operators. For this reasons 
we shall not perform such a study in our model-independent analysis. A study in a model- 
dependent context considering both SM and heavy gauge bosons contributing to the B 
X s j branching ratio is in progress [45]. 



3 New contributions at the high scale 
3.1 Effective Hamiltonian for b — > S7 

Adopting the overall normalisation of the SM effective Hamiltonian and considering as a 
first step only the dipole operators and the contribution of new physics to their Wilson 
coefficients, we write the effective Hamiltonian relevant for b — > 57 at the high matching 
scale hh 3> Mw as: 



,(6-^7) _ 4G F l 
V2 



H ef{ ~ 7z V ts V tb 



AC 7 ^fx H )Q 7 ^+AC 8G (fM H )Q 8 G+^C' 7j (^ H )Q' 7j +AC' 8G (fj, H )Q' 8G 

(3.1) 

At this scale, the W bosons are still dynamical and the usual SM contribution to the 
Wilson coefficients is absent. We will include it after integrating out the W boson at the 
lower scale ~ 0(i%). In our conventions the dipole operators are given by: 

l07T z 

(3.2) 

Qsg = ^2^^ ^ Pr bp G% 

and the primed operators are obtained from them after replacing the right-handed projec- 
tor, Pr, with the left-handed one, Pl- In the SM the contributions of the primed operators 
are suppressed by ra s jra\, relative to those coming from Q 77 and Qsg- We decompose the 
Wilson coefficients at the scale nw as the sum of the SM contribution and the new one from 
the exchange of neutral flavour gauge bosons, after evolving the new physics contribution 
using the RG running of the operators Qip and Q^q, namely 

Ci{fiw) = Cf M ( m ) + AC^w) (3.3) 
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and similarly for the primed coefficients. The SM Wilson coefficients at LO are [46] 



nSM (li \ _ 3x f ~ 2x t i 8xf + 5xf - 7x t 

°^ {m) ~ 4(x t - l) 4 lnXi " 24(x t - 1)3 ' (3 - 4) 

Csg 0*0 - 4(xt _ 1)4 " 8(xt _ 1)3 . ( 3 - 5 ) 

where xt = rrif/My^. The corresponding primed coefficients are also given by the previous 
expressions, but with an extra suppression factor of m s /mb. 

The new physics contributions at jiw^ ACi(fj,w), follow from the RG evolution of the 
Wilson coefficients at hh, ACi(fj,n), down to fiw', Their relations are discussed in Sect. 4. 
We can further decompose the new physics contributions at the high scale as the sum of 
two pieces, one coming from the exchanges of heavy exotic quarks and the second from the 
exchanges of light SM down-quarks: respectively, 

AC 7j (m) = AC 7 h 7 eavy (^) + AC^ h \v H ) , 

(3.6) 

ACsgM = AC^{ m ) + ACf(^) . 
In the following we analyse these two contributions separately. 

3.2 Contributions from heavy exotic quark exchanges 

We begin with the new contributions to b — > sj arising from the exchange of heavy neutral 
flavour gauge bosons and heavy exotic quarks. In what follows, we restrict the discussion 
to the contribution from a single virtual gauge boson An and a single virtual heavy quark 
of electric charge —1/3 and mass m' F . The results will then be generalised to an arbitrary 
number of such particles. Prom the general Feynman rules listed in Eqs. (2.14) and (2.15), 
we read off the flavour changing couplings of Ah to the bottom and strange quarks and 
the exotic quark F relevant for b — > sj: 



< = i 9 -fl» (PlC[ F + P R C f / 

(3.7) 

The only diagram 1 with a virtual F and Ajj exchange contributing to the b — > S7 transition 
in the unitary gauge is shown in Fig. 1. In the same figure, we also show the analogous 
diagram contributing to b — > s g. The latter will also contribute to the B — > Xgj rate via 
QCD mixing of Q%g i n Qi^- in both diagrams, there is no triple gauge boson vertex as the 
heavy gauge boson Ah considered does not carry electric or colour charge. 




1 The corrections on the external legs, necessary for a canonical kinetic term, are understood to be 
contained in finite off-diagonal field renormalisations. 
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1,9 



Figure 1. Magnetic penguin diagram contributing to b — ¥ sj and b — > sg with internal Ah and F 
exchanges. 

We further decompose the Wilson coefficients AC 4 heavy at the high scale as the sum 
of the SM-like LL contribution and a new LR one, where L (R) stands for the Pl (Pr) 
projector in the vertex of Eq. (3.7) involving the s(6)-quark: 

ACf7( m ) = A LL C 7 h 7 cavy (/^) + A LR C^( m ) , 

(3.8) 

AC 8 h r y (M = A LL C^(u, H ) + A LR C^( m ) . 

The corresponding primed coefficients to Eq. (3.8) are obtained by interchanging L -H- R 
and including the suppression factor m s /mb in the primed Wilson coefficients. 

Following [35], the SM-like contributions LL and RR are obtained by noting that their 
topological structure is the same as the SM gluon magnetic penguin. Consequently, the 
known SM results for the Wilson coefficient of the gluon dipole operator, C^q , can be 
used to obtain both LL and RR contributions after a suitable redefinition of couplings 
and masses. Notice that in the neutral exchange there is generally no analogue of a GIM 
mechanism, therefore mass independent terms in the SM Wilson coefficients must be kept. 
We extract them from [47, 48]. Our results are as follows: 

ALL^heavy/ s _ 9}J 1V1 W U L °L / r< SM (^ _L 1 

^ ^--S£W H ^SV^\ CaG{X) + V' (3 . 9) 



with x = m? F /M\ H . The RR contributions are obtained by the interchange L f-> R. 
C^q (x) is the SM function on the right-hand side of Eq. (3.5). 

We now consider the LR and RL contributions, which have no SM equivalent. To 
obtain the Wilson coefficients, we adapt to our case known calculations of b — > sg in the 
context of b — > 57 in the left-right symmetric models [34, 47], where, however, the process 
is mediated by charged gauge bosons. In this way, we obtain the following LR Wilson 
coefficients: 

1 f ,2 M 2 I fisF* r bF 

A ^ M --ZtiW H ^^^ * G{h (3.IO) 

A M C heavy (Mff) = ^LR^yy^ _ 

with 
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The analogous RL contributions are obtained by interchanging L -R- R. In this case 



The following properties should be noted: 

1) the LR (RL) contribution in Eq. (3.10) dominates over the LL (RR) one of Eq. (3.9) 
in a large part of the parameter space, due to the factor m' F /mb; 

2) the factor m! F /rrib is replaced by mt/m\, in the usual left-right symmetric models, 
therefore in these models the enhancement is less pronounced; 

3) C^q{x) is a non- vanishing monotonic function of x and takes values in the range 
[—1, —1/4] for x from to oo; 

4) in the decoupling limit these contributions turn out to vanish. Indeed when large 
masses for the exotic fermions fn'p k are considered, the masses of the neutral gauge 
bosons Ma h are also large, as we have already commented in Sect. 2. Therefore the 
LL (RR) contributions approach zero with 1/M\ and the LR (RL) ones are strongly 
suppressed by m' F /M\ H . 

3.3 Contributions from light SM quark exchanges 

We now discuss AC^ ht and ACgg ht from Eq. (3.6), namely the contributions arising from 
the exchange of light down- type quarks in loop diagrams analogous to the one in Fig 1, 
where the flavour changing couplings of Ah to the SM quarks read 



arising from the general Feynman rules listed in Eq. (2.12). In the SM, the light-quark 
contributions cancel each other at scales well above their masses due to the GIM mechanism. 
This is in contrast to the present case, where such a mechanism is absent. 

Following the procedure of the previous section, we decompose the Wilson coefficients 
at the high scale as the sum of a LL part and a LR part: 



A H 




(3.12) 



(3.13) 



In the framework of effective field theories the light degrees of freedom are treated as 
massless at the high matching-scale \xh- Therefore, we only need to account for the first 
term in the expansion of the light masses and only mass-independent terms in Eq. (3.9) 



and (3.10) contribute. For the LL SM-like contributions we find 

ALL^light/ x _ 1 9h 1V1 W V^°L °L f 1 
A W 7 W - "ft 172- 2^ 



e 5 2 2 Mi w ^ y; s y ifc y^/ (314) 



and for the LR contribution we have from Eqs. (3.10) and (3.11) 



77 M ~ M| Mi H ^6 ( 1} ' (315) 



where / stands for the SM down-type quarks. Analogously, we obtain also the results for 
the primed Wilson coefficients. The LL and RR SM-like contributions have already been 
considered in [49] , while the LR and RL are new. Notice that the latter come with a factor 
of mf/mb and hence, the effects from d- and s-quarks are suppressed with respect to the 
6-quark contribution. 

3.4 Generalisations and the GIM-like mechanism 

It is easy to extend the previous results to the case of an arbitrary number of neutral flavour 
gauge bosons, Ar., and heavy fermions, F^, by performing the following substitution in 
Eqs. (3.9) and (3.10): 



9h n sF* ribF ( r<SM(„\ i 1\ , \ ~* ^Hj n ski* n bki ( r<SM(„ \ , ^ 

M 2 U L,R U L,R y^8G { x ) "•" g J > U L,R \^8G \ x ki) T g 



m'. 



(3.16) 



9h 



2 " b F r< sF * r< hF n L R(^ v \ " „ 2 _Ej> r<ski* n bki r>LR(„ \ 

M 2 U L,R U R,L °8G \ x > > Z^ 9H ^ lf2 U L,R °i?,,L °8G \ x ki) > 



where x^i = m' 2 p k /Mj, and are the couplings among the light quarks (s, b) with 

if j J 

Aff i and F^. In this substitution we take Ajj i to be the mass eigenstates of the neutral 
gauge bosons. Similarly, we may also generalise Eq. (3.14) and (3.15). 

Notice that in the case that the masses of the heavy particles span a wide range of 
values, the generalisation above is only a first approximation: a more precise result follows 
from a detailed analysis using the results of Sect. 4 considering several threshold scales. 
This further study goes beyond the scope of the present paper and in the following we only 
consider a single high matching scale. 

The factor m' F /M\ can represent a dangerous enhancement of the LR (RL) con- 
tribution in Eq. (3.16), which is unlikely to vanish, barring accidental cancellations. We 
identify a series of conditions under which this contribution exactly vanishes: 

1) 9Hi = 9H-, i-e. there is only one gauge symmetry or all the new gauge coupling 
constants have the same strength; 
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2) m' Fk = m' F , Ma h _ = Ma h , i.e. all the exotic fermions and the neutral gauge bosons 
are degenerate in mass, as in the case of unbroken flavour symmetries; 

3) Yli C s lr Cth, = J2i ^ s lr Cff l = 0, i.e. the couplings are unitary matrices satisfying 
Cl = Cr. 

Similar conditions apply also to the contributions of SM down- type quarks. 

These conditions assure the exact cancellation of the NP contributions. However, 
similarly to the GIM mechanism, if they are only partially satisfied, the terms in Eq. (3.16) 
do not cancel each other completely and the NP contributions could still be dangerous. 

Furthermore, we remark that condition 2) would correspond to degenerate SM fermions 
and thus these conditions cannot be satisfied in models which successfully explain the SM 
fermion mass pattern through the see-saw mechanism illustrated in Sect. 2. In these models, 
however, the mixings between SM and exotic fermions are small, cf. Eqs. (2.14) and (2.15). 
This provides a sufficiently strong suppression to safely neglect these contributions. We 
further discuss this point in Sect. 5.2. 

4 QCD corrections 
4.1 General structure 

In order to complete the analysis of the b — > s'j decay we have to include QCD corrections, 
which play a very important role within the SM, enhancing the rate by a factor of 2 — 3 [50]. 
It originates dominantly from the mixing of charged current-current operators into the 
dipole operators and to a smaller extent from the mixing with QCD-penguin operators. 

When the flavour gauge bosons and the exotic fermions are integrated out at 
the high scale fj,H, in addition to the dipole operators discussed at length in the previous 
section, also neutral current-current operators corresponding to a tree-level Ajj t exchange 
are generated together with ordinary QCD- and QED-penguin operators. 

The contributions of QCD- and QED-penguin operators, arising from diagrams with 
An i and exchanges, are much less important than within the SM. Indeed, in the SM 
context, QCD-penguins do not have to compete with tree- level FCNC diagrams, as the 
latter are absent due to the GIM mechanism. However, in the present case there are 
neutral current-current operators originating from the tree-level Au i exchanges that are 
not suppressed by a s (fj,n) and loop effects, contrary to QCD-penguins. Thus, for all 
practical purposes the contributions of QCD-penguins from Ajj i and exchanges at the 
high scale \xn can be neglected. Even less important are effects from QED-penguins. Note 
that this is in contrast to dipole operators which cannot be generated at tree-level, but 
nevertheless mediate the b — > s"f decay. 

In what follows we will extend the SM RG analysis by considering the QCD effects of 
neutral current-current operators at scales fif, < fj, < fin generated by diagrams with the 
exchange of a single neutral gauge boson Ah- The extension to the case of an arbitrary 
number of such gauge bosons is straightforward. This analysis does take into account 
the mixing under QCD renormalisation of these neutral current-current operators 1) into 
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dipole operators, 2) among themselves, 3) into QCD-penguin operators. The latter mixing 
generates contributions to QCD-penguin operators at scales lower than lih, even if the 
initial conditions at lih of the Wilson coefficients of these operators are neglected. Even 
if this effect is numerically negligible, we include it for completeness as the mixing of 
QCD-penguin operators and dipole operators is taken into account in the SM part. 

Before going into details, let us note that the neutral current-current operators do not 
mix with the charged current-current operators as their flavour structures differ from each 
other. On the other hand, similarly to charged current-current operators, neutral current- 
current operators have an impact on dipole operators and QCD-penguins through RG 
effects without being themselves affected by the presence of these two types of operators. 

We would like to emphasise that the QCD analysis presented below is also relevant for 
other processes, such as non-leptonic two-body B decays, and other observables, like e'/e. 

Denoting symbolically charged current-current, QCD-penguin, dipole, and neutral 
current-current operators and the corresponding primed operators respectively by 



Q cc i Qp , Qd , Q nn , Q'p , Q'd j Q T 

the structure of the one-loop anomalous dimension matrix looks as follows: 



(4.1) 





Q cc 


Qp 


Qd 


Q nn 


Q cc 




x 2 


X 3 





Qp 





x A 


x 5 





Qd 








x 6 





Q nn 





Yx 


Y 2 


^3 



(4.2) 





x 4 


x 5 





Q'p 







x 6 





Q'd 




Yy 


Y 2 


Y 3 


Qnn I 



Q'p Q'd Q nn ' 

The non- vanishing entries denoted by Xi are known from the SM analysis [51, 52]. The 
new results in the present paper are the entries denoted by Yi, the initial conditions for 
Qd, given already in Sect. 3, and the corresponding conditions for Q nn , given at the end 
of this section. 

The RG analysis of the NP contributions can be performed independent from the SM 
one: that is 



Cl rp (ji b ) = W ij (ji b , m )Cf p ( m ) 



(4.3) 



where Uij and Wij 



are the elements of the RG evolution matrix. The Wilson coefficients for 
the primed operators at the high scale are in general different from the expression above, 
but the evolution matrices U and W do not change. The general structure of these matrices 
in terms of the anomalous dimension matrices and the QCD-/3-functions is well known [53] 
and will not be repeated here. Therefore, we present in the next sections our results for 
the matrices Yi and the initial conditions for Q nn and Q nnl operators. 
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4.2 Operator basis for neutral current-current operators 

There are 48 neutral current-current operators generated by the exchange of Ah, all con- 
taining the neutral currents (s 7^ Pl,r b) and the flavour conserving currents (/ 7^ Pl,r f) 
with f = u,c, t, d, s, b. At this scale also the exotic fermions have been already integrated 

OUt at fJLH- 

The general notation for the 48 Q nn and Q nnl operators in question will be 

Q f 12 (A,B), for A, B = {L, R} . (4.4) 

I^or 1 hs*Ij stucc 

Q\{L, R) = (s a 7m P L bp) (up ^ P R u a ) , 

(4.5) 

Ql(L, R) = (s a 7m P L b a )(up 7" P R up) , 

where a, /3 = 1, 2, 3 are colour indices. 

Moreover there are 8 additional neutral current-current operators that have to be in- 
cluded, as they have a different structure, when colour and flavour structures are considered 
simultaneously: 

Qi(A, B) = (s a llx P A dp) (dp r Pb b a ) , 

M) 

Qt(A,B) = (s alfl P A d a )(dp^P B bp) . 

In this classification, we refer to Q nn [Q nn '] as those operators Qi^(A,B) with A = L 
[A = R] and B = L, R [B = R, L). 

We point out that the operator basis in Eqs. (4.4) and (4.6) can be reduced using 
Fierz transformations. Moreover some relations can be found between the usual charged 
current-current and QCD-penguin operators 

Qi = (s a 7m Pl cp) (cp 7 M P L b a ) , 
Q2 = (s a 7 M Pl c a ) (cp 7 M P L bp) , 

Qz = (s a 7m Pl b a ) Yl fe> ^ Pl q P) ' 

q=u,d,s,c,b 

Qa = (s a 7m Pl bp) ]T (qp ^ P L q a ) , (4.7) 

q=u,d,s,c,b 

Qb = (s a 7m p l b a ) ^2 {qp 7 M Pr qp) , 

q=u,d,s,c,b 

Qa = (s a 7m p L bp) Yl ^ Pr *0 • 

q=u,d,s,c,b 

Yet, similarly to [54, 55] we have decided to work with all operators and not use such 
relations when performing the RG analysis in order to keep the anomalous dimensions in a 
transparent form. As discussed in [54, 55] one can work with linearly dependent operators 
in the process of RG analysis without any problems and use such relations only at the end 
of the evolution if necessary. However, in the case of CV 7 (//&) at LO this is not required. 
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In this context we should warn the reader that at NLO level a use of Fierz relations in 
the RG analysis could lead to wrong results within the NDR scheme unless Fierz evanescent 
operators are introduced [56, 57]. 

4.3 Anomalous dimension matrices 

We define the LO anomalous dimension matrix through 

TK) = ^7 (0) , (4-8) 

and give below the results for 7^ (dropping the index (0) to simplify the notation). 

The inspection of the one-loop diagrams contributing to the anomalous dimension 
matrices shows that only 16 operators have to be considered in order to find the full 
matrix. These are 

Ql 2 (L,L), Q( 2 (L,L) , Qi 2 (L,L), Q d lj2 (L,L), 
Qi 2 (L,R), Q l2 (L,R), Q s l2 (L,R), Q 12 (L,R). 

The anomalous dimensions for u replaced by c and t are equal to the one of Q\ 2 . The 
ones of Q\ 2 are equal to the ones of Q\ 2 . The anomalous dimensions of the remaining 28 
operators, namely the primed operators obtained by L f-> R, are the same as those of the 
corresponding unprimed operators (see Eq. (4.2)). 

The mixing between the Q nn operators and Qr>, that is the matrix Y 2 , can be extracted 
from [51, 58] by inspecting the mixing between QCD-penguin and dipole operators. For 
the transposed matrices we get 





Qi 


Ql 


Qi 


Qi 


Ql 


Q s 2 


Qi 


Qi 




416 
~8f 





232 





232 
~~8T 


232 
~"8f 





232 
~~81 


Q&G 


70 
27 


3 


70 
27 


3 


151 

~W 


151 

"27~ 


3 


70 
27 



(4.10) 



1d(L,R) 





Qi 


Ql 


Qi 


Qi 


Ql 


Q2 


Qi 


Qi 


Ql-y 


448 





200 





200 


16 


80 


32 


~~8f 


"81~ 


"81~ 


~9 


~Y 


~~9 


Q8G 


119 


-3 


119 


-3 


173 


16 


-4 


8 












3 



(4.11) 



The mixing between neutral current-current operators, that is the matrix I3, is uni- 
versally given by two 2x2 matrices: 



T n (L,L) 





Qi 


Q2 






Qi 


Q2 


Qi 


-2 


6 


j nn (L,R) = 


Qi 


-16 





Q2 


6 


-2 




Q2 


-6 


2 



(4.12) 
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Only operators with the same flavour-conserving structure (/ 7^ Pl,r f) mix in this way. 

Finally, in the case of Y\ there is a universal mixing of Q nn operators into QCD-penguin 
operators: 







Q3 


Q 4 


Q5 


Qe 


IP - 




2 

~9 


2 
3 


2 

~9 


2 
3 



(4.13) 



for 

Q nn = ^ Q u,c,t,d,s,b^ Q s,b^ ^ Qd^ ^ Q u,c,t,d,s,b ^ R ^ ^ 

with no mixing for the remaining operators in Eq. (4.9). 
4.4 Initial conditions 

The initial conditions at hh for the dipole operators have already been presented in Sect. 3: 
see Eqs. (3.4), (3.5), (3.9), (3.10), (3.14) and (3.15). Here we give the corresponding initial 
conditions for the neutral current-current operators. 

The initial conditions are determined by integrating out all the heavy degrees of free- 
dom at fin and matching the result to the effective theory. At LO this matching is trivial. 
In the following we give the result for the simplified case in which only one neutral flavour 
gauge boson Ah contributes to the initial conditions. The generalisation to several gauge 
bosons is obvious. 

Denoting the vertex f{Ajjfj as in Eqs. (2.11) and (2.12), the general expression for 
the Wilson coefficient of the neutral current-current operator 

Q((A, B) = (s a 7ai P A b a )(f p ^ P B f p ) , (4.15) 

in the normalisation of Eq. (3.1), is given by 

A AB C^ H ) = 

while all the coefficients A AB C\ are zero 
operator is 

A AB c d 2 M - 

while A AB Cf (fin ) is again zero. 

5 Phenomenological analysis 
5.1 The B — > X s j branching ratio 

The SM prediction for the B — > X s 7 branching ratio at NNLO [50, 59, 60] reads, 

Br(B -> X s j) = (3.15 ± 0.23) x 10~ 4 , (5.1) 



1 „2 tut2 (~<sb* /-iff 

.2 A#2 

Iff 



(4.16) 



2giMi H V t * s V tb 
Similarly, the Wilson coefficient for the Q%(A, B) 

1 9h M 2 W Cg C b B d 
2g\M\ H V t * s V tb ' 
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and has been calculated for a photon-energy cut-off E~ > 1.6 GeV in the £?-meson rest 
frame. This is to be compared with the current experimental value [61], 

Br(B X s7 ) = (3.55 ± 0.24 ± 0.09) x 10~ 4 , (5.2) 

for the same energy cut-off E~. In the presence of NP the expression for the branching 
ratio is given as follows: 

Br(B -> X sl ) = R (|C 77 (/i b )| 2 + |C^> 6 )| 2 + N{EJ) , (5.3) 

where R = 2.47 x 10~ 3 and N(E 1 ) = (3.6 ± 0.6) x 10~ 3 . The parameter R is simply 
an overall factor, whose determination is discussed in [59, 60], while N(E^) is a non- 
perturbative contribution. Calculating the NP contributions in the LO approximation, 
but including the SM contribution at the NNLO level we find for m = 2.5 GeV 

C 7j (fM b ) = C^ M (fi b ) + ACV 7 ( M6 ) (5.4) 

where the central value of CjJ^ 1 ' (nb) at the NNLO, corresponding to the value in Eq. (5.1), 
is given by [50, 59, 60] 

Cffifib) = -0.3523 (5.5) 
and the NP one, calculated by us at the LO, by 

AC^ifib) = K7 ACV 7 0u ff ) + k 8 AC 8G (/xh)+ 

+ £ k{ a A^ClOM + E^A A LA Cf fa). (5-6) 

A=L,R A=L,R 
f=u,c,t,d,s,b 

Here k's are the NP magic numbers listed in Tab. 1, calculated taking a s {Mz = 91.1876 GeV) = 
0.118. The primed coefficient C' 7ry (nb) can be obtained from (5.4)-(5.6), by interchanging 
L <-> R and taking the initial conditions of the primed Wilson coefficients. In particular, 
the NP magic numbers listed in Tab. 1 are also valid for the primed case, as QCD is blind 
to the fermion chirality 2 . 

The SM contribution containing NLO and NNLO QCD corrections exhibits a negligi- 
ble /^-dependence. This is not the case for the NP contribution at LO. However, we have 
checked numerically that when the NP contribution enhances the SM value by 20%, the 
//^-dependence in the total branching ratio amounts to a 3% uncertainty for /i& S [2.5, 5] 
GeV. For smaller deviations from the SM prediction the uncertainty further decreases. 
This renders this uncertainty sufficiently small for our purposes. 

Before concluding a few observations are in order: 

2 In the computation of G^(jj,b) we are neglecting the SM contributions Cf^" J '(fj,b) and Cgcf ' (pb), and 
the NP contributions AC7 7 (/X{,) and ACgail^b), since all of them are suppressed by m a /m,b with respect 
the unprimed ones. Our numerical analyses confirm that this approximation is consistent with errors below 
the percent level. 
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m 


200 GeV 


1 TeV 


5 TeV 


10 TeV 


M z 


K 7 
K 8 


0.524 
0.118 


0.457 
0.125 


0.408 
0.129 


0.390 
0.130 


0.566 
0.111 


u,c 
K LL 

K LL 

LL 

8,6 

K LL 

"d 
LL 


0.039 
-0.002 
-0.040 
0.087 
0.128 


0.057 
-0.003 
-0.057 
0.090 
0.147 


0.076 
-0.002 
-0.072 
0.090 
0.163 


0.084 
-0.001 
-0.079 
0.090 
0.168 


0.030 

-0.032 
0.084 
0.116 


u,c 
K LR 

K LR 
K d 

LR. 

s,b 
K LR 
~d 

LR 


0.085 
0.004 
-0.015 
-0.078 
0.473 


0.128 
0.012 
-0.025 
-0.092 
0.665 


0.173 
0.023 
-0.036 
-0.106 
0.865 


0.193 
0.028 
-0.041 
-0.111 
0.953 


0.065 

-0.011 
-0.070 
0.383 



Table 1. The NP magic numbers for ACV 7 (//f,) defined in Eq. (5.6). For completeness, we include 
in the last column the case of a flavour- violating Z. The analogous table for ACsaiHb) is provided 
in Appendix A. 

1) Similarly to the SM, the magic numbers K7 and Kg suppress the initial values ACV 7 (//#) 
and ACsg(^h)- This is due to the QCD RG evolution running down to /ib- Further- 
more, the suppression of AC , 7 7 ( / u//) increases with ^h- 

2) As in the SM, provided A AE 'C^ifJ-u) and A ab C2(^h) are sufficiently larger than 
AC7 7 (//#), the additive QCD corrections stemming from the mixing of the neutral 
current-current operators into the dipole operators are dominant. Furthermore the 
QCD factors increase in most cases with ^h- The most prominent is the coefficient 
k d LR which could even be of 0(1), but also k^ r , k d LL and k s £ l lr are sizable. However, 
values of the corresponding initial conditions could compensate these coefficients, as 
in the case of small couplings among SM fermions. 

3) The SM and NP primed dipole Wilson coefficients, C[{nw) and AC'^h)-, are sup- 
pressed by m s /mb and turn out to be numerically negligible. On the other hand, 
this suppression is absent in the contributions of the primed neutral current-current 
operators Q nnl and therefore they should be considered in the determination of the 
branching ratio. 

5.2 Three classes of models 

The results for Kj shown in Tab. 1 are model-independent and hold for all models in which 
the neutral gauge bosons have flavour violating couplings to fermions as discussed in Sect. 2. 
On the other hand, the initial conditions of the Wilson coefficients are model-dependent. 
However, in spite of the large variety of models it is possible to distinguish three main 
classes, even if hybrid situations are also possible: 
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Models without exotic fermions. This is the case of models that do not aim at 
explaining the SM fermion masses and mixings or models that do provide such an 
explanation but without exotic fermions. For instance, Z' constructions in which the 
theory is anomaly free without extending the fermion spectrum of the SM fall into 
this class. This is possible when the Green-Schwarz mechanism is implemented in 
the theory or when the generator of the additional U(l) is a linear combination of 
the hypercharge Y and B — L, where B is the baryon number and L is the lepton 
number. 

For all models of this class, the new diagrams with virtual exotic fermions do not 
contribute. On the other hand, the contributions from the exchange of the light 
down-type quarks, the neutral current-current operators and the corresponding QCD 
evolution are still present. As illustrated in our numerical analysis of Sect. 5.4, the 
light-quark contributions turn out to be negligible in the whole parameter space. 

Models in which the SM fermion-mass patterns are governed by exotic fermions 
through a see-saw mechanism as illustrated in Sect. 2. In this case, the couplings of 
the neutral gauge bosons to SM and exotic fermions, described by Eqs. (2.12)-(2.15), 
are suppressed by sin0L h ,R k given in Eqs. (2.7) and (2.8). In the specific limit of 
m' ft 3> rrif k and ~ = M D the suppression in sin^ fc) /j fc is approximately of 

order O (yJm fk /m' F ^ . 

In order to illustrate the size of the contributions from the new diagrams of Sect. 3, 
we redefine the couplings C s ^ l R and C^ R to exhibit the dependence on the suppressing 
factor. Without loss of generality, we can approximately write 




Ci k R^J^Ci%, Cf^J^Cfy (5.7) 



nip 



As a result, for the case of arbitrary Ajj t and F^, the expressions on the right-hand 
side of Eqs. (3.9) and (3.10) are given by 



^ W) ""6?^M1- rn> V t * s V tb { C °° M + 3 



A 



(5.8) 



Lfi^hcavy/ s + \ ^ m W m s °L ° R nLR(„ 



and A( LL < LR ^Csc Vy (u, H ) = -3A( LL < LR ^C^ avy {u, H )- Similar expressions can be writ- 
ten for the primed contributions. Notice that the m' F /mi enhancement is completely 
removed from the LR (RL) Wilson coefficient and is replaced by the suppressing fac- 
tor ^rris/mb. The only dependence on m' Fk is in the loop factor C^Q^Xi^). On the 
other hand, an extra inverse power of m'p k is now present in the LL (RR) contribution, 
which further suppresses this term. 

As we shall explicitly demonstrate in Sect. 5.4 these contributions turn out to be 
negligible with respect to those from the neutral current-current operators. The 
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same holds for the light-quark contributions. Hence, when dealing with this class of 
models, we simply neglect the contributions of exotic fermions and therefore, as in 
1), only neutral current-current operators and their mixing with dipole operators are 
relevant. 

3) Models with exotic fermions, in which the definition of their couplings to gauge bosons 
and SM fermions is independent of the mechanism of the SM fermion mass generation. 
For example, this is the case of models in which these couplings do not arise from the 
standard kinetic terms of Eq. (2.9). As a result, sm0L k ^ k can in general be much 
larger than in the previous classes of models. 

In this case, the expressions in Eqs. (3.9) and (3.10) do not suffer from the additional 
suppressions of Eq. (5.8) and in particular the LR (RL) contribution is strongly 
enhanced by the factor m' F /rrn,. For the models of this class, the contribution of 
exotic fermions is the dominant one and all the other contributions can be safely 
neglected. 



We summarise the relevant features of the models in these three classes in Tab. 2. 



Classes of models 




qT 


Q nn , Q nn ' 


1) without exotic fermions 

2) with exotic fermions and see-saw 

3) with exotic fermions but without see-saw 


Absent 
Negligible 
Dominant 


Negligible 
Negligible 
Negligible 


Dominant 
Dominant 
Negligible 



Table 2. Summary of the different classes of models and the corresponding NP contributions to 
b — > sj: contributions from the exchange of heavy exotic quarks Q h ^ avy , from the exchange of SM 
down-type quarks Q l ^ ht ', and from neutral current- current operators Q nn and Q nnl '. 



5.3 Model-independent constraints 

It has been pointed out in [62-67] that considering the measurements of B — )■ X s l + £~ 
[68, 69] the sign of CV 7 (//&) in the presence of NP is likely to be the same as in the SM. 
This provides a first model- independent constraint for NP contributions: since C^(nb) 
is negative, ACy 7 (/ib) should also be negative in order to soften the tension between the 
central values of the SM prediction and the experimental determination of Br(B — > Xgj). 
This holds not only for contributions from unprimed operators, but also when considering 
the primed ones due to the relative suppression of C' 7 (fj,b) with respect to CV 7 (/ib), as 
explicitly confirmed by our numerical analysis. 

The sign of the NP contributions is determined by the product of the initial conditions 
and the QCD magic numbers in Tab. 1. The signs of the K{ factors are fixed solely by the 
QCD running, while the sign of the initial conditions ACj(^#) depends on the couplings 
of the new neutral gauge bosons to the fermions. 

In particular, for a model in the first two classes, this constraint translates into the 
requirement that the second line on the right-hand side of Eq. (5.6) should be negative, 
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namely that: 

AC 77 (/i fe )^ Y, 4a^ LA cIM+ y £k d LA A LA C$(n H )<0, (5.9) 

A=L,R A=L,R 
f=u,c,t,d,s,b 

where only the couplings listed in Eq. (2.12) are involved. 

On the other hand, when a model belongs to the third class, the first line on the 
right-hand side of Eq. (5.6) must be negative, 

AC 77 (yu 6 ) ~ k 7 AC 7i (hh) + k 8 AC 8G ((m h ) < . (5.10) 

This puts a constraint on the couplings of Eqs. (2.14) and (2.15). 

Once the coupling constants are chosen such that the NP contributions have the same 
sign as the SM one, it is possible to further constrain the parameter space by requiring that 
the predicted branching ratio should not exceed the experimental bound. From Eq. (5.3) 
we find the constraint 

- AC 77 (/x b ) + 1.4 (|AC 77 (/i b )| 2 + |AC 77 Gu 6 )| 2 ) < 4.2(6.1) x 1(T 2 , (5.11) 

corresponding to the lex (2 a) departure from the experimental value, Eq. (5.2). 

It is straightforward to apply these constraints to all models belonging to one of the 
three classes discussed above. It is not possible to obtain more insight without specifying a 
particular model (an analysis will follow in [45]). However, in the next section, we provide 
a simplified representative for each class of models and show how the constraints discussed 
so far apply to each of them. 

5.4 Toy-model examples 

For each class discussed above, we consider a toy-model in order to justify the approx- 
imations made in the previous sections and to exemplify the application of the model- 
independent constraints. 

Classes 1) and 2): Here, the relevant initial conditions are those presented in Eqs. (4.16) 
and (4.17). The coupling constants entering these expressions are C^ R , Cff R , Cff R 
and C^ R , with f = u,c, t, d, s, b. To simplify the analysis, we assume that all flavour- 
violating couplings with a strange flavour are equal to C S FV , the two bottom couplings 
Cj^ R equal to C FV and all the flavour conserving ones equal to Cpc- Furthermore, 
we take Vt s = —0.04047, Vtb = 0.999146 [70], assume that gn = 92-, and consider only 
one heavy neutral gauge boson. In this way we have defined a toy model with only 
four free parameters: three coupling constants and the mass of the neutral gauge 
boson. 

In Fig. 2 on the left we show the breakdown of C 77 (/ib) in its different contributions as 
a function of the coupling constants Cfv = C FV = C FV = Cfc, for Ma h = 1 TeV. 
For completeness, we also show the exotic-quark contributions for m' F = 10 TeV. 
As expected from the discussion in Sec. 5.2, all the NP contributions apart from the 
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neutral current-current are negligible, justifying our approximations in the previous 
section. We remark that Cy (//;,) almost coincides with the neutral current-current 
contribution (red line). In Fig. 2 on the right we show the value of the coupling 
constants CpVi as a function of Ma h , for which the bound in Eq. (5.11) is saturated 
at the la and 2a level. As we can see, for small values of Ma h the couplings are 
constrained to small values. 




0.2 0.4 0.6 0.8 1 200 300 500 1000 2000 

\Cfv\ M AH [GeV] 

Figure 2. On the left, the different contributions to CV 7 (/i&) (solid line) are plotted as functions 
of Cfv'- i n blue the SM contribution, in red the neutral current- current contribution, in green the 
overlapping contributions from exotic and SM down-type quarks. The shadowed region is excluded 
imposing the bound in Eq. (5.11) at la (lighter) and 2a (darker). On the right, we show the value 
of Cfv as o, function of Ma h for which the bound in Eq. (5.11) is saturated. Again the shadowed 
region represents the excluded values. 

In Fig. 3 we separately present the implementation of the two model-independent 
constraints, for Ma h = 1 TeV. On the left, the constraint on the sign of ACj^dif,) 
mostly reduces the parameter space to cases in which the coupling constants C S FV 
and Cp V have opposite signs. In Fig. 3 on the right, the bound from the branching 
ratio also applies, further constraining the parameter space. 

Class 3): We consider here the case with only one heavy neutral gauge boson and one 
exotic fermion. The relevant initial conditions are those in Eqs. (3.9) and (3.10). 
The coupling constants which enter these expressions are Cff R and C^^. Fixing 
9h = 92, Ma h = 1 TeV and m' F = 10 TeV, and identifying the coupling constants 
Cfv = C'jfft = Cjf R , we illustrate in Fig. 4 on the left that the dominant NP con- 
tributions stems solely from exotic quarks. For a complete comparison, we also plot 
the contributions from light quarks and neutral current-current operators, adopting 
the same conventions as for Fig. 2. 

In Fig. 4 on the right we show the value for the coupling constants Cfv as function 
of Ma h for which the la and 2a bounds in Eq. (5.11) are saturated. As we can see 
the constraint on the couplings is very strong and only a small part of the shown 
parameter space survives even for large values of Ma h - 
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Figure 3. On the left (right), the points satisfying ACV 7 (/it;,) < (Eq. (5.11)) in the plane C FV 
vs. Cp V (= Cpc). The shadowed regions now represent the values for which the bounds are passed: 
lighter (darker) intensity refers to —4.2(6.1) x 10~ 4 . 



0.5 
0.4 



3 °' 3 

O 0.2 

I 




0.1 

F~ 








0.2 



0.4 0.6 
\C FV \ 



0.8 



- 



0.8 



0.6 



0.4 



0.2 








200 300 500 1000 
M AH [GeV] 



2000 



Figure 4. A similar description to Fig. 2 applies. Here, however, the green (purple) line refers to 
the exotic-quark (SM-quark) contributions. 



In Fig. 5 on the left (right) we show how the couplings are constrained by the ex- 
pression in Eq. (5.10) (Eq. (5.11)): a negative value for AC7 7 ( / u; ) ) is only recovered 
when both Cff and C b j^ R have the same sign, while the lower bound in Eq. (5.11) 
provides a very strong constraint on the parameters of the model. 



6 Conclusions 



Extensions of the SM in which the gauge group is enlarged by additional symmetries 
are attractive, because they provide an explanation of the flavour puzzle and predictions 
testable at colliders. In these models heavy exotic fermions are usually introduced in order 
to cancel possible anomalies and to justify the observed SM fermion spectrum through the 
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Figure 5. On the left (right), the points satisfying ACV 7 (/if,) < (Eq. (5.11)) in the plane Cff vs. 
CW , for C h jf € [—1, 1]. The shadowed regions are defined in Fig. 3. 



see-saw mechanism. On the other hand, the presence of new heavy neutral gauge bosons 
and exotic fermions in principle translates into NP contributions to FCNC processes. 

In this paper, we have pointed out two new contributions to the B —¥ X s ^y decay that 
arise in such models. The first is generated through one- loop diagrams mediated by neutral 
gauge bosons and exotic down-type quarks. The relevance of this contribution depends on 
the strength of the left-handed and right-handed flavour violating couplings of the neutral 
gauge bosons to the SM and the exotic quarks. Analogous effects can be present in other 
processes like [/, — > ej and t — > cy as well as flavour conserving observables like (g — 2)^ 
and EDMs, in which dipole operators play the dominant role. 

The second contribution is due to the presence of neutral current-current operators, 
mediated by neutral gauge bosons, and arises only if the neutral gauge bosons have flavour 
violating couplings to the SM quarks. Through the QCD mixing with the magnetic dipole 
operator Q7 7 , these new neutral current-current operators contribute to B — > X s j. To our 
knowledge these QCD effects have been calculated here for the first time. Furthermore, 
our QCD analysis of the mixing among neutral current-current, QCD penguin and gluonic 
dipole operators could also be relevant for other processes, such as non-leptonic two-body 
B decays, and other observables, like e'/e. 

Beside these new contributions, we have also considered the contributions arising from 
one-loop diagrams with the exchange of neutral gauge bosons and SM down-type quarks, 
that have been only partially analysed in the literature. 

We have studied the impact of all these contributions in a model-independent approach 
and summarised the resulting constraints in Eqs. (5.9)-(5.11). We have presented these 
expressions in such a manner that in order to test a specific model, it is sufficient to specify 
only the values of the couplings of the neutral gauge bosons to SM and exotic quarks and 
their masses. In particular, it is not necessary to repeat the QCD analysis. 

A detailed application of this study on a concrete NP scenario is in progress. Here, 
without entering into details of a particular model, we have described three representative 
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classes of models and discussed the relevance of these contributions. For models in the 
first class, the SM spectrum is enriched only by the neutral gauge bosons, but no exotic 
quarks are present. In this case the contributions from the exotic quarks are absent, but 
those from the neutral current-current operators turn out to have a potentially observable 
effect on the branching ratio of B — > X s j. The value of the masses of the neutral gauge 
bosons and the strength of their flavour violating couplings to SM fermions determine the 
relevance of this effect. 

The second class accounts for models in which the SM fermion masses m / are explained 
through the see-saw mechanism with heavy exotic fermions of masses m' F . In this case, 
the couplings of neutral gauge bosons to SM and exotic fermions are suppressed by terms 
mf/m' F . Therefore, the contributions to Br{B —> Xg'y) from the exchange of exotic quarks 
turn out to be negligible with respect to those from the QCD mixing of the neutral current- 
current and magnetic dipole operators. 

The models in the third class are characterised by the presence of heavy exotic fermions, 
which either provide the SM fermion masses through a different mechanism than the see- 
saw or do not participate at all in the explanation of the SM flavour spectrum. In this case, 
no suppression occurs in the couplings of neutral gauge bosons to SM and exotic fermions 
and the contributions to Br(B —> X s ~f) from exotic quarks are enhanced by terms mfp/rrib, 
where nib is the mass of the bottom quark. These contributions dominate over all the 
others. 

For all the models in the three classes, the contributions from SM down-type quarks 
turn out to be subdominant in the whole parameter space. 

Our analysis shows once more how FCNC processes can put constraints on beyond- SM 
constructions even before the discovery of new particles in high-energy processes. 
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A The B — »■ X s g decay 

Similarly to Eq. (5.4) we also evolve ACsg down to fib = 2.5 GeV to obtain 



ACgjifib) 



p 7 ACV 7 (/i ff ) + p 8 AC 8G (/iff)+ 



+ 53 P [ A a la c((vl h ) + Y.pla & LA ciM 



(A.l) 



A=L,R A=L,R 
f=u,c,t,d,s,b 



with the NP magic numbers pi listed in Tab. 3. 
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m 


200 Gev 


1 lev 


5 lev 


10 lev 


M z 


P7 

















P8 


0.568 


0.504 


0.456 


0.439 


0.607 


u,c 

Pll 


-0.124 


-0.138 


-0.147 


-0.150 


-0.115 


Pll 


-0.015 


-0.033 


-0.046 


-0.050 


_ 


Pll 


-0.124 


-0.138 


-0.147 


-0.150 


-0.115 


s,b 

Pll 


-0.243 


-0.279 


-0.307 


-0.318 


-0.222 


PL 


-0.119 


-0.141 


-0.160 


-0.168 


-0.107 


u,c 

Plr 


0.184 


0.229 


0.270 


0.287 


0.160 


Plr 


0.015 


0.037 


0.055 


0.062 




Plr 


0.184 


0.229 


0.270 


0.287 


0.160 


8,6 

Plr 


0.311 


0.382 


0.447 


0.474 


0.272 


Plr 


-0.064 


-0.052 


-0.034 


-0.025 


-0.067 



Table 3. The NP magic numbers for ACscil^b) in Eq. (A.l). For completeness, we include in the 
last column the case of a flavour-violating Z. 
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